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Abstract 

We study the propagation of time-harmonic acoustic or transverse magnetic 
(TM) polarized electromagnetic waves in a periodic waveguide lying in the half- 
strip (0,oo) x (0,L). It is shown that there exists a Riesz basis of the space of 
solutions to the time-harmonic wave equation such that the translation opera- 
tor shifting a function by one periodicity length to the left is represented by an 
infinite Jordan matrix which contains at most a finite number of Jordan blocks 
of size > 1. Moreover, the Dirichlet-, Neumann- and mixed traces of this Riesz 
basis on the left boundary also form a Riesz basis. Both the cases of frequen- 
cies in a band gap and frequencies in the spectrum and a variety of boundary 
conditions on the top and bottom are considered. 
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1. Introduction 

Periodic media have received much attention in recent years since they can 
prohibit the propagation of electromagnetic and acoustic waves in some fre- 
quency ranges [l[ . This localization property is a consequence of band structure 
of the spectrum of the underlying differential operator and of the presence of 
band gaps in this spectrum. Waves with frequencies in a band gap will decrease 
exponentially inside such media as a consequence of the exponential decay of 
the Green's kernel, seej2L The band structure of the spectrum is explained 
by Floquct theory, see [Jill, [H[ . We refer the reader to p for a mathematical 
introduction to photonic crystals and to [7, 8] for a physical introduction. 
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The strong localization property can be used to construct devices that mould 
the flow of light or sound at a very small length scale. The simulation of such de- 
vices requires the numerical solution of differential equations in locally perturbed 
periodic media, which is a challenging task. The proof of existence and numer- 
ical computation of defect modes has been studied in 0, E3, EH, 0, 0, 13, 14 1. 
Problems in locally perturbed infinite periodic media with a source term consid- 
ered as part of the problem have been studied more recently in [IB , T^, 17 1 . Here 
a natural approach consists in solving a boundary value problem on a compact 
set enclosing the perturbation and using a Dirichlet-to-Neumann or a related 
operator on the artificial boundary. Hence, the problem is decomposed into an 
interior and an exterior problem. The purpose of this paper is to contribute to 
the understanding of exterior boundary value problems in semi-infinite periodic 
waveguides. A summary of our results has already been given in the abstract, 
for a precise formulation we refer to the next section. 

The study of wave propagation in doubly periodic half planes can be re- 
duced via the Floquct transform to wave propagation in semi-infinite periodic 
waveguides with quasi-periodic boundary conditions on the lateral boundaries. 
Moreover, using a clever trick proposed by Fliss & Joly [TBI, boundary value 
problems in the doubly periodic exterior of a square can be reduced in some 
sense to boundary value problems in a doubly-periodic half plane. Therefore, 
the analysis of this paper is also relevant for these problems. 

Our results prove two open conjectures formulated in connection with a 
numerical method to compute the Neumann-to-Dirichlet map proposed by Fliss, 
Joly and Li (see the discussion of Theorem l2.2[) . Actually, this study arose from 
attempts to justify an alternative numerical approach, which will be published 
elsewhere. Moreover, our results explain with a new approach the exponential 
decay of waves in periodic media with frequency in the band gap studied in 
Tol . HI 0, EH, EH and even provides the optimal decay rate which corresponds 



to the one of the slowest decaying Floquet mode, which may provide guidance 
for photonic crystal optimization by focusing on this Floquet mode and trying 
to make its decay as fast as possible. 

The paper is outlined as follows: in § [3] we introduce the problem and present 
the main theorem. Some known prerequisites are collected in § [3J in particular a 
generalization of Rouche's theorem shown in 18( , which will serve as an essential 
tool in the following analysis. The remaining part of the paper is dedicated to 
the proof of the main theorem and will be summarized at the end of section [3j 
There are two appendices on radiation conditions and on uniqueness results. 



2. Statement of the problem and the main results 

The propagation of time-harmonic acoustic or transverse magnetic (TM) 
polarized electromagnetic waves in a 2-D waveguide lying in the half-strip S + := 
M + x (0, L) is described by the differential equation 

Av + lo 2 s p v = mS+. (la) 
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We assume that e p <E L°°(S) with S:=lx (0, L) is periodic with period length 
1 in the first variable and bounded away from 0, i.e. 

£ p (xi + 1,2:2) = £ P (xi,x 2 ), for all (xi,x 2 ) € S, 
< essinfep < e p < e a.e. in 5* 

for some e > 0. On the top and bottom of S + wc consider a boundary condition 
L 
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Figure 1: Dielectric permittivity of the semi-infinite waveguide. 

yv{x 1 ,-) = (0,0) T , x x > (lb) 

with one of the following boundary value operators 7 : H 2 ((0,L)) — > R 2 : 

Dirichlct -y u v := (v(0), v(L)) T 

Neumann := (v'(0), v'(L)) T 

mixed 7dn^ := (i>(0), u'(L)) T 

/3-quasi-periodic 7,91; := (e l/3 u(0) - w(L), e f V(0) - «'(L)) T , /3 € [0, 2tt). 

To treat /3-quasi-periodic boundary conditions with /3 = and /3 — n we have 
to impose the symmetry condition e p (xi,X2) = e p (xi, L — X2), x € (0, 1) x (0, L) 
for reasons explained in § [6j 

Furthermore, we imposed a boundary condition on the left: 

dv 

rv = f with tv :=e B v(0 r ) + 9 N — (0,-). (lc) 

ax 1 

Here #d,#n S C with |6*d| + |(?n| > 0, and we consider r as an operator with 
values in a space which depends on #n and 7 and will be defined in § 

To describe our condition on the behavior of the solution as x± — > 00, which 
will complete the formulation of the boundary value problem, we need the fol- 
lowing definition: 

Definition 2.1. A Floquet mode is a nontrivial solution to (Tla)) , (| lb|) of the 

form 

m 

exp(i£xi)22x™~ j Uj(xi,X2) 
3=0 

with £ 6 C and functions itj satisfying + 1,2:2) = ^(2:1,2:2) for all 

(2:1,2:2) e 5+. We call £ the quasi momentum and m the order of the Flo- 
quet mode, assuming that u^ m ^ ^ 0. 
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Note that the quasi momentum is only determined up to an additive integer 
multiple of 27r. Typically we will choose € [— 7r,7r). We complete the formu- 
lation of the boundary value problem by the radiation condition 

v G H^ + (S + ) where H} f - + (S + ) := 1^(S + ) © span{v+, • • • ,«+}, (Id) 

where wj 1 ", . . . , are Floquet modes with real quasi momentum, the choice of 
which is described in the following: 

It is known (see [l9|, Corollary 5.1.5] or Corollary 14. 6j) that the wave guide 
supports a finite even number 2n of linearly independent Floquet modes with 
real quasi momentum. To formulate a radiation condition we introduce the 
scsquilincar form 



q Xl (v,w) := J (-^-(x)w(x)-v(x)^-(x)jdx2, xi>0 (2) 

for solutions v,w to (|la,|) and (|lb[) . Since q Xl is actually independent of x\ 
as a consequence of Green's theorem, we omit the index x\ in the following. 
If the time dependence is given by exp(— iujt) with uj > 0, then Qq(v,v) is 
proportional (with positive constant) to the energy flux through a cross section 
{xi} x (0,L) (see e.g. [3, § 5.6.3] where the sign convention exp(icji) is used). 
From a physical mode we expect that energy is transported to the right, so 
Qq(v,v) > 0. It can be shown (see [lj| Theorem 5.3.2]) that there exists a 
basis {v^ , • • • , i>i, • • • , v^} of the span of all Floquet modes with real quasi 
momentum, which consists of Floquet modes satisfying the orthogonality and 
normalization conditions 

q( v j~> v k) = > for j,fce (3a) 

l( v t> v k) = i6 J,k, q(vj ,Ufe) = -i6j,k, for j, fee {l,...,n}.(3b) 

(As discussed in Appendix |Appcndix A[ spanji;^, ■ ■ ■ ,v£} is not necessarily 
uniquely determined by this condition for all uj, but the following results hold 
for any choice of the 

Each element v € H} ( ' + (S + ) has a unique representation of the form v = 
v + 2™=i a nV„ with v € H} / (S + ) and a n <E C, and we introduce a norm on 

by ||»||^, + (s+) = \\v\\ 2 myi s +) + ELi Kl 2 - 
Let us introduce the translation operator 

(Tv)(x) := v(x\ + 1, a; 2 ) 

and assume that the Floquet modes v+ have been chosen such that H} ( ' + (S + ) is 
invariant under T- (This is always the case if uf, . . . , «i are all of order 0.) Let 
C -ffy + (S + ) be the linear subspace of all weak solutions to (|la|) and (|lb[) . If 
the boundary value problem ((T|) is well posed, i.e. the operator r|y : V — > fll 
has a bounded inverse (see Prop. 16.21 and Appendix |Appendix B| for sufficient 
conditions), we can define the monodromy operator 

TZ := rT(rly)- 1 : -> fl}, 
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which maps the trace of a solution to ([1} to its trace at one periodicity length 
to the right. We are now in a position to formulate our main result: 

Theorem 2.2. 1. There exist Floquet modes v£ , n S N the first n of which 
are defined as described above, which form a Riesz basis of V . This basis 
can be chosen such that T : V — > V is represented by an infinite Jordan 
matrix, which contains at most a finite number of Jordan blocks of size 
greater than 1 all of which are of finite size. 
2. If problem ([1]) is well posed for some choice of t, then {tv+ : n € N} 
is a Riesz basis of H'L, and with respect to this basis the operator 1Z is 
represented by the same Jordan matrix as T in part\^ 

The second part of this theorem proves the conjectures in [ItI Remark 5.1] 
and [3 Conjecture 3.2.52]. (I n the latter case, to prove that ^| S p an ^ TU +. n >^ 
has spectral radius < 1, we also have to take into account Proposition 14.31 ) 

Wc point out that the wcll-poscdncss assumption in the second part of Theo- 
rcm l2.2l is always satisfied for certain types of Robin trace operators r (see Prop. 
16.21 and IB.2| . We will use this fact in the proof of the first part of the theorem. 
For other trace operators r non-uniqueness may occur at certain frequencies ui. 



3. Preliminaries and outline of the proof 

3.1. Sobolev spaces on the cross section 

For the boundary value operators 7 <E {7d, 7n, 7dn, 7^} defined in the intro- 
duction, we define the second derivative operators 

L> 7 : V(D-y) -> i 2 ((0, L)), w i-> -w" 

with V{D 1 ) := {w e H 2 {(0,L)) : 7 w = (0,0) T }. 

It is well known that these operators are positive and self-adjoint with compact 
resolvents. Moreover, complete orthonormal systems of eigcnpairs 
{(^fc 7 ^> (^i 7 "*) 2 ) '■ k *eT} are known explicitly: 

boundary condition (t) k£' T 

Dirichlet 4/f sin(^i) ^ N 

Neumann ^\ cos (?£t) ^ NU{0} 

/3-quasi-periodic yj± exp (j£±™it) Z 

To simplify our notation we will often suppress the dependence of and itj. 7 ^ 
on 7 in the following. Sobolev spaces corresponding a boundary value operator 
7 G {7d,7n,7dn,7/3} can then be defined by fl"'((0,L)) = V((I + D 7 ) s / 2 ) 
with norm : = II (I + D^) s ' 2 w\\l,2 for s > 0, and for s < the space 



mixed J isin ll±^lt Z^pz N 
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H*((0,L)) can be denned as completion of L 2 ((0,L)) under the norm || • \\h»- 
More explicitly, 

\\wf H , = ^2(i + kly\(w,4> k )f (4) 

kei 

It will be convenient to renumber the square roots of the k;'s in increasing 
order such that 

a(D 7 ) = {k\ : k G 2} = {rf : I G N} 

and < Ki < K2 < . . • • This defines a bijective mapping N — > I, k i— > £"(£;). 
Note that 

271 \ , 27r , 

K 2 z+i = Ki + — / and K2/+2 = «2 + -f-l (5) 

for 1 = 1,.... Moreover, we will write tpk — ipuk) f° r k G N. 

5.12. Sobolev spaces on the strip 

We introduce the self-adjoint negative Laplace operators 

L> 7 : £>(L> 7 ) -> L 2 (S), u^-Aw, 

with D(I> 7 ) := {w G H 2 {S) : Van G R7«(i 1 ,-) = (0,0) T } 

for 7 G {7d, 7n, 7dn, 7^}- Then we define 

^(5):=2?((/ + J D 7 ) s / 2 ), S >0, 

and := ||(1 + D 7 ) s ^ 2 v\\ L 2. For s < 0, H*(S) is defined as completion 

of L 2 (S) under the norm || ■ \\h"- It can be shown that 11^(8) = Hq(S), 
H^(S) = H^S), and H^(S) = {v G H^S) : e^v(-,0) = v(-,L)}, and the 
norm || • \\ H i is equivalent to || • \\ H i given by ||v||^i = J s (|w| 2 + \S7v\ 2 )dx. 
H^(S + ) can be defined as the set of all restrictions to S + of functions in H^(S) 
with ||w||_ff s (s+) := m f{||^lli?4(S) : v\s+ = v}. Moreover, the trace operators 

r D :H^S + )^H^ 2 ((0,L)), v^v(0,-), 

r N : H} f (S+; A) H^ 2 ((0, L)), v ^ ^(0, ■), 

are well defined, continuous, and surjective, and have bounded right-inverses 
(see (H). (Here IMI^ (S+ . A) := \\v\\ 2 H i (s+) + || A v || 2 2(5+) .) We choose H; := 

H] /2 ((0,L)) if N =0 and RZ < := H~ 1/2 ((0, L)) if 6> N ^ 0. 

3.3. Floquet transform 

The Floquet transform is defined by 

T : L 2 (S) -> L 2 ((-7r,^),L 2 (fi)) 



.Fu(a, a;) := V v( Xl + I, x 2 )e- ia[xi 
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where f2 = R/Z x (0,L). It is isometric and its inverse is given by v(x) = 
-^J^Tv(a,x)e^ da, x € S (@). 

We will frequently use the orthonormal bases {tpm,n : m £ Z,n £ N} of 
L 2 (fi) defined by 

:= exp(27rimxi)V' i ( 2 7) (a;2), x £ f2 . (6) 

If w 6 Hi(S), then for all a € [— 7r,7r] the function J r v(a,-) belongs to the 
Sobolev space defined by H*(f2) := {u £ L 2 (£l) : ||u||#»(n) < 00} with 

IMk(n) : =( E (l + l«| 2 )1(^ (7) ,«)| 2 ) 1/2 

for s > and as completion of L 2 (Q) under this norm for s < (@). For a G 
[— 7r,7r] the operator a£ 7 ' : H 2 (fl) —> L 2 (Q) uniquely defined by the property 

A a (Fv(a,-)) = {FAv){a,-) 

is given cxplicitely by 

A Q = e - MXl Ae^ 1 = (<9 X1 + iaf + <9 2 2 = A + 2ia<9 xi - a 2 . (7) 

5.^. Floquet modes and characteristic values of (B%) 

Due to ((7J the mapping a i— > A Q is a polynomial with coefficients in 
£(H s (il), iJ s_2 (Sl)) for all s, so in particular, it has a holomorphic extension 
denoted by 

A e := A + 2i£8 Xl -i 2 , (eC. (8) 
Moreover, let us introduce the operators 

S e : -> H°- 2 (n) B^u := A s u + cu 2 e p u 

for C, which are well defined for any s £ [0,2]. For a review of the properties 
of holomorphic extensions of Floquet transformed periodic differential operators 
in a much greater generality we refer to 

£0 £ C is called a characteristic value of £ 1— > £?£ if is not injective. If 
£0 is a characteristic value of B$ for some choice of s £ [0,2], and B^ Uo = 
for ito £ H*(£l) \ {0}, then uq £ H 2 (£l) by elliptic regularity results, and £0 is 
a characteristic value of B^ for any parameter s £ [0,2]. Therefore, the set of 
characteristic values does not depend on the parameter s, and for studying this 
set we may choose s at our convenience. The operators B^ are defined such that 
(with s = 2): 

Remark 3.1. The following statements are equivalent for £0 £ C: 

• £0 is a characteristic value of (B^). 

• There exists a Floquet mode v(x) = e^ 0X1 u(x) with u £ H 2 (il). 
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It follows from the second statement that: 

Remark 3.2. If £ g C is a characteristic value of £ H> B^, then £ + 27rZ is a 
characteristic value for all Z € Z. 

To formulate some results by Gohberg and Sigal [l8| in Theorem 13.41 below, 
which will serve as an essential tool in the following analysis, we first have 
to recall the definition of the multiplicity of characteristic values. Note that 
for all £ the operator is a compact perturbation of the operator A — 1 in 
C(H^(Q), L 2 (f2)), and that A — 1 has a bounded inverse. Therefore, by analytic 
Fredholm theory [2lJ the set of characteristic values of B^ is discrete. For the 
special case at hand the definitions in [l8[ simplify as follows: 

Definition 3.3. Let 6i and 62 be Banach spaces and £ 1— > B^ a holomorphic 
mapping defined on a domain in C with values in £(61,62). Moreover, let 
B^ = B + where B has a bounded inverse in £(62,61) and K% is compact 
for all £. 

The point £0 G C is called a characteristic value of (5g) if there exists a 
holomorphic function £ i— > (called root function) with values in 61 such that 
it£ 7^ and B^ Q uc a = 0. (Note that such root functions, which may be chosen 
constant, exist if and only if B^ a is not injective.) The multiplicity of a root 
function (u^) is the order of £0 as a root of £ n- B^u^. u G 61 is called an 
eigenvector of (6^) corresponding to £0 if u = uj for some root function (u^) 
of (B^) corresponding to £o- The ranA; rank(u) of an eigenvector u is defined 
as the maximum of all the multiplicities of root functions (u^) with u = U£ . 
(Under the given assumptions the geometric multiplicity a := dimker6^ of £0 
and the ranks of all eigenvectors u € ker6^ are finite, see 0, Lemma 2.1].) 

A canonical system of eigenvectors of (6^) corresponding to £0 is defined as 
a basis {w 1 ' ,...,u( a >} of ker6j with the following properties: rank?/ 1 ) is the 
maximum of the ranks of all eigenvectors corresponding to £01 and ranku^ 
for j = 2, ...,a is the maximum of the ranks of all eigenvectors in some direct 
complement of spanjw/ 1 ), u^ -1 '} in kcr6^ . The numbers Tj = ranku^ 
(J = 1, a) are called the partial null multiplicities of the characteristic value 

and n((6^); £0) = ( r i, r 2, r a ) the a-tuple of partial null multiplicities. We 
call 9T((-Bf ); £0) = ?*i+r2 + ...+r a the (total) null multiplicity of the characteristic 
value £0 01 (-£>£)• 

Finally, we call canonical Jordan chains associated to the characteristic 
value £0 an y sets {uq , ii£ _i}j 1 < i < a, of vectors in 61 such that 
{uq 1 \ Wq"^} is a canonical system of eigenvectors and £0 is a root of order 
r s of £ ^ ^(EfcLo 1 ^ ~ Co)^) for any 1 < j < a. 

Let r be a simple, closed, rectifiable contour contained in the domain of 
analyticity of (65) and of (B^)^ 1 and let £1, £ n be the characteristic values 
of (6^) enclosed by T. (Recall that simple means that T has a continuous, 
bijective paramctrization 7 : [0, w) — > T.) Then we set 

n 

9t((6 c );r):=5>((6 e );^). (9) 

3=1 
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We cite the following generalization of Rouche's Theorem: 

Theorem 3.4 (flijl). Assume that (A%) satisfies the assumptions of Definition 
\3.3\ and let T be a simple, closed, rectifiable contour in the domain of analyticity 
of (Aj) and of (A^ 1 . If £ h-> is another holomorphic function defined on 
the same domain as {A^) with values in L(B\,B2) and if 

\Al x St&<\ for all £er, 

then (Aj + ) _1 is analytic in some neighborhood ofT and 

m((A ( + S ( );T)=m((A i );T). (10) 

We can now state the following generalization of Remark 13.11 



Proposition 3.5. Suppose that £o is a characteristic value of (B^) with par- 
tial null multiplicities n((-Bf);£o) = (r\,r2, ■■■,r a ). Then, the vector space Vj 
of all Floquet modes with quasi momentum £o is a direct sum of T -invariant 
subspaces V^ j of dimensions rj G N, j — 1, ...,a. Moreover, there exists a 
basis {vjfi, . . . ,Vj tr .-i} ofV^ _j with respect to which T : V^ a j — > V^ .j is repre- 
sented by a single rj x rj Jordan block, and Vj^ are Floquet modes of order k. 
Moreover, in the expansion 

Vj , k (x) = e*°* ]T ITZW"^)' k = 0,..., rj -l (11) 
i=o [ >■ 

we have Uq' ^ = ■■■ = M^'-i ^ = > an d ■ ■ ■ ? ^o*^} * s a basis °f 

ker(B 5o ). 

Proof. From [lj| § 3.4.3] we know that there exist canonical Jordan chains 
{uq i , G*.' 7 '_ 1 }, 1 < j < a associated to £o, and Vf = ©" =1 V£ ,j where V^ 0J = 
spanju^fc : < k < rj — 1} and 



Since 



(ix 1 +i) k l _ Uh 
(k-l)l 

1 i m {ix 1 ) k - l ~ m 



^ ^ to! (k - I - to)! 1 

1=0 m=0 v ' 

--e^±^e^ k fp^^u^(x) 
^ to! ^ (k-m-l)\ ' 



jk—m 



Z *° ^ (fc-TO)!^'' m(a;) ' 



m=0 
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the subspaces Vf 0) j are T-invariant and T is represented by an upper triangular 
Toeplitz matrix M^ j with respect to the basis (uj,fc)o<fc<rj-i of V^ j. Since 
the entries of the first upper diagonal of M^ a j do not vanish, we have (M^ j — 
e ifo/y = if and only if I > rj, so M^ j is similar to a Jordan block Jj of size 
rj with diagonal values exp(i£o)- It is easy to see by induction in Tj that this 
similarity transform can be achieved by an upper triangular matrix D with l's 
on the diagonal, i.e. D~ M^ jD = Jj. Then Vj^ '■— Vj,k + z2i=o ^l,k^j,l is a 
Floquet mode of order k of the form (jlip with u^'^ = Uq f , and T is represented 
by Jj on the subspace Vf 0) j with respect to the basis {'Uj Oj ■ ■ • ? w j,r,-i} of V^ 0l3 -. 

□ 

5.5. Outline of the proof of Theorem ] 2. Oft 

We amend the finite number of right propagating Floquet modes v^, . . . , ti 
which have been described in section[5]by an infinite number of decaying Floquet 
modes v+, n = n+ l,n + 2,..., which are chosen according to Proposition 13.51 
for each characteristic value of (Bg) with positive imaginary part. These i?+ 
are arranged in increasing order of the imaginary part of the corresponding 
characteristic values and will be properly normalized. The set {v£ : n £ N} is 
a Riesz basis of V if and only if the operator 

oo 

T : Z 2 (N) -> H} f -+(S + ) 7 (a n ) ^ £ a„v+ (13) 

n=l 

is a norm isomorphism from Z 2 (N) to V. Our strategy is to compare T with a 
reference operator T r corresponding to the case uj = 0, which we will refer to 
as the unperturbed case. In the unperturbed case all interesting quantities can 
easily be computed analytically. We will proceed as follows: 

• In section 2] we establish the existence of a countable number of character- 
istic values and derive precise estimates on the difference of characteristic 
values with large imaginary parts in the perturbed and the unperturbed 
case using Theorem 13.41 

• In scction[5]wc show that T—T T is compact by estimating the perturbation 
of eigenvectors of (B%) for £ with large imaginary parts. Moreover, we show 
that T r is a norm isomorphism from l 2 (N) to its range. 

• In scction[B]we show for trace operators r, which satisfy the well-posedness 
assumption in Theorem l2.2l that rT r : l 2 (N) — > H'L is a norm isomorphism. 
Moreover, we show injectivity of T. Together with Ricsz theory and the 
well-posedness assumption this readily implies that both tT : I 2 (N) -4 
and T : I (N — > V are norm isomorphisms. Moreover, by construction (see 
Proposition 13. 5|) both the operators TZ and T are represented by Jordan 
matrices. 
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4. Estimates of characteristic values 



4-1. Characteristic values for u> — 

Recall from £13.11 that K n are the eigenvalues of the negative Laplacian on 
the cross section with boundary conditions 7 and set 

£m,n : = — 27TTO + lK n , 771 6 Z, tlfN. 

Lemma 4.1. The characteristic values o/(A^) are precisely the numbers ^ m ,n 
andt; m ^ n (counted with their total multiplicities if these numbers are not distinct) 
with m G Z and n € N. If K n > aZZ partial null multiplicities are 1 and i/ 
K n = t/ie partial null multiplicity is 2. ^4n eigenvector corresponding to both 
£,m,n and £, m ,n is given by the function (p m ^ n defined in ((6]). 

Proof. Note that —Aip m! „ = ((27rm) 2 + K$)<p m ,n = \€m,n\ 2 <Pm,n for (m,n) € 
Z x N and define u(m, n) := (u, (f m ,n)- Due to ([8]) we have 

-(A?T)(I) = -(A + 2z£(2™m) - £ 2 )^(0 = (|6| 2 + 2^m + £ 2 )£(Z) 
= (£-£l)(£-5)«(0> i=(m,n)eZxN. 

Since {ipi : I E Z x N} is an orthonormal basis of L 2 (f2), this shows that the 
problem separates into a countable set of scalar problems indexed by I = (m, n), 
each of which has the two simple characteristic values and if n n > and 
the characteristic value of multiplicity 2 if K n = 0. □ 

Since the characteristic values are always 27r-periodic, also for ui > (see 
Remark l3.2j) . it suffices to study the characteristic values in the strip {z € C : 
— 7r < 5ftz < 7r}. Therefore we define 

£j = ik„ = £ ,n, i£N 

as "reference characteristic values". Note that £Jj corresponds to the physi- 
cal exponentially decaying Floquet mode exp(i^xi) = exp(— K n x\) if n n > 
whereas £^ corresponds to the unphysical exponentially growing Floquet mode 
exp(^xi) = exp(K n xi). 

4-2. Estimates of "large" characteristic values for Dirichlet boundary conditions 

For the sake of clarity we first prove estimates on the location of character- 
istic values for the case of Dirichlet boundary conditions 'y — 'yrj where K n — ^p- 
before treating general boundary conditions. The main tool will be Theorem 

] 

We define open discs V n and closed rectangles lZ n by 

V n :=!z€C:\z-a< - ^ j, (15) 
Tin ■= [z 6 C : < 7T,\%Z-K n \ < . (16) 
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The radii of these discs are chosen such that || A^ 1 uj 2 e p \\c(L 2 (n)) < 1 for au 
£ £ dV n as shown later. Moreover, we choose N <E N such that all disks T> n 
with n > N are mutually disjoint and contained in the strip {z € C : |5Rz| < 7r}. 
i.e. 

f k„ + cj 2 emax(l, L) 1 , , 

TV > min < n £ N : 2 > I ' ^ ^ 

Proposition 4.2. For Dirichlet boundary conditions we have 

m((B^),dD n ) = 1 = m({B^,dR n ) for alln>N. 

In particular, all characteristic values of (B^) in the half-strip {z 6 C : — n < 
< 7r, 5z > kat — 7r/(2X)} are simple and contained in one of the disks T> n . 

Proof. For each n > N we apply Theorem[33]with = A 6 e £(£ 2 (0), #~ 2 (0)) 
and Sc = B^ — = cj 2 Ep. We will show that 

II VeplU^p,)) < 1 £ e ft„ \£>„ . (18) 

Then Theorem 13.41 is applicable with T = &D n and yields 0T((-Bg), 9Z?„) = 
^((Aj), 9D„) = 1 (with the help of Lemma T4.ip . and it follows from the Neu- 
mann series or another application of Theorem [33] with T — d!Z n that there are 
no characteristic values of (B$) in 7Z„ \ T> n . This proves all claims. 
Since \\w 2 e p \\c(L 2 (n)) = u 2 s, eq. (fT8j) holds true if we can show that 

5 we 
By virtue of (fT4|) this is equivalent to 

supl 1 :leZxN|<4--. £eft„\2V 

ll£-6ll£-6l J w£ 

Note that 

TZ n = {£ e C : VZ e Z x N |£ - &,«| < |£ - 61} , (19) 

so |C - 61 > If - 6j,«| > Kn 2 fnLi for a11 £ e ^» \ T °g ether witn thc 
inequality |f — 61 > l^f — 961 > ^£ > |(k„ + K n -i) we obtain the desired 
estimate 

|£ - 61 If - 61 >^ 2 £ for all I 6 Z x N 



(see Fig. 2(a) I. Since the supremum over I above is obviously attained, we get 
a strict inequality. □ 
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4-3. Estimates of "large" characteristic values for general boundary conditions 
For general lateral boundary conditions, in particular quasi-periodic bound- 
ary conditions, the eigenvalues K n may come in pairs which are arbitrarily close 
together or coincide. Recall, however, that we have K n +2 = K n + j for all 
n £ N and for all trace operators 7. We assign to each index n its group G(n) 
as follows: G(n) := {n} if K n +i — n n = -j-, i.e. if the K n are equidistant, which 
is the case for 7 £ {7d,7n,7dn}- If i^n+i — K n < |j, then G{n) := {n,n + 1} 
and if n n+ i — K n > j-, then G(n) := {n,max(n — 1, 1)}. 
To each group we assign a set 

7J G(n) :={^C:V/eZxN min |£ - £| < |£ - , (20) 

n'eG(n) 



which is a closed rectangle since the form a rectangular grid (see Fig. 2(b) ). As 
remarked in (fT9)) , these rectangles coincide with the ones defined in the previous 
subsection if 7 = 7rj. Moreover, we introduce open disks, which again coincide 
with those of the previous subsection for 7 = 7d: 



V n := {zeC:\z-C\ < 7 



u 2 e 



Proposition 4.3. Define the set 
S := I J V n with iV > min/n £ N : inf > - £max ( 1 ' L ) \ 

(21) 

('see [£(&2p. Then all characteristic values of (-Bf) m {z £ C : SRz £ 
[— 7r, 7r], Sz > inf2 g i{ G{JV) 3z} are contained in S and for each connected compo- 
nent S of S we have 

dS) = <tt((A € ); &S) = #{n £ N : £ £ S} . (22) 

Proof. The proof is analogous to that of Proposition ^. 21 Note that N is chosen 
such that S C {z £ C : |3z| < 7r} and UneG(n) ^« c ^G(n) f° r au n > N. 
By the same arguments we can show that 1 1 1 oj 2 e p ||/;(i 2 (o)) < 1 f° r all £ £ 
ft G(ri) \SifScU„ eG( „)A,. □ 

4-4- Number of "small" characteristic values 

Only for very small uj the argument of Proposition l4.3l can work for all n since 
for large u> some of the discs T>„ overlap with their neighbors. In the following 
we will determine the number of characteristic values with "small" imaginary 
parts. First we cite the following symmetry result for the set of characteristic 
values: 

Proposition 4.4 (0, Theorem 5.3]). For real-valued e the set of character- 
istic values of (Be) is symmetric with respect to the real axis, and moreover 
m((B s ); Co) = !Tt((B € ); £0) for all £ € C. 
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-2vr 



2n 



(a) Proof of Proposition 14, 21 



(b) Proposition [33] 



Figure 2: Panel (a): An important ingredient of the proof dictating our definition of the 
geometry are uniform lower bounds on products of the distances indicated by fat double 
arrows. 

Panel (b): The shaded area shows the set S defined in (1 2 1 I t in the case of quasi-periodic 
boundary conditions. The crosses indicate the characteristic values of (B^), and the dots 
indicate some of the reference characteristic values (A^). The dashed lines show the boundaries 
of rectangles TZG(n) defined in H20I I. 



To count the number of these "small" characteristic values in the strip {z £ 
C : — 7r < 5R(z) < 7r} we surround them by a contour containing none of them. 
By Proposition 14.31 the segment [zn, zjy + 2ir] with zjy := — ^ + * i n fzeTC G(JV) ^ z 
contains no characteristic value. 

Proposition 4.5. Choose zjv as above and any complex path P from zn to zn 
such that r :— [zn + 27r, zn] Li PL) [z~n, zn + 2tt] Li — (2ir + P) encloses precisely all 
characteristic values of (B^) in the rectanglelZ := {z £ C : 5Rz £ [—w,tv), \5sz\ < 



\Qzn\} (see Fig. 3(a)). Then 

n((B s y,r) = 2N. 

Proof. Let us define -B(£,/i) := + fie p for /i £ [0,w 2 ]. From Proposition! 
we can deduce that the segments [zn',zn + 2it] and [z/v + 2n, zjv] contain no 
characteristic values of £ H> B(£,fi) for any fi £ [0,w 2 ]. For a given [i £ [0,u; 2 ] 
it follows from the discreteness of the set of characteristic values of £ h-> B(£, /j.) 
that the segment [zjv,ztv] contains at most a finite number of characteristic 
values, and there exists S > such that 5-balls around these characteristic 
values contain no further characteristic values. We deform the straight line 
[zn,Zn] hi (^-neighborhoods of the characteristic values to left semi-circles as 



shown in Fig. 3(a) to obtain a contour P^. Because of the 2n periodicity of 
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Figure 3: Panel (a): The dashed line indicated the integration contour in the proof of Propo- 
sition [43J Crosses indicate characteristic values of (B^), circles mark four of the "reference" 
characteristic values of (Aj). 

Panel (b) : The proof of Lemma 15,21 involves lower bounds of the distances of points £ in a 
connected component S of the set S (the dark shaded region) to all "reference" characteristic 
values £m,n not contained in S, see 



the characteristic values, the imaginary parts of the characteristic values on 
[zn, ~zn] and on 2tt + [zn, zjv] coincide, and replace the segment 2ir + \zn, zn] 
by — (2tt + Pfj). This yields a contour T M := [zn + 27r, zn] U P m U [zn, zm + 2n] U 
— (27r + Pfi) which contains precisely all characteristic values of B(-,n) in the 
rectangle 7Z. Moreover, T w 2 coincides with the contour T in the proposition. 
With this construction the function N : [0,u; 2 ] — > {0, 1, 2, . . . }, 

N(p) :=«n((fl f );IV) 

is well-defined, and in particular it does not depend on the choice of . More- 
over, it follows from Theorem 13.41 with A% = B(£,fi) and S$ = 5e p such that 
\S\ < 1/ (e max^ e r t , || j 4^ 1 ||£(L 2 (si)) ) that N is constant in a neighborhood of 
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each fi £ [0,u! 2 ]. Hence, N is constant on whole interval [0,w 2 ] and 
3t((fl e ); T) - N(u 2 ) = N(0) = Ot((A f ); T ) = 27V. 

□ 

We mention that Proposition l4.5l vields an independent proof of the following 
well-known result (see [ljj Corollary 5.1.5]): 

Corollary 4.6. The dimension of the space V p spanned by Floquet modes with 
real quasi momentum is finite and even. 

Proof. Due to Proposition l4.4l the contour T in Proposition l4.5l cncloses the same 
number rn of characteristic values with positive and with negative imaginary 
parts. Together with Proposition 13.51 we find that dimV p = 2N — 2m is finite 
and even. □ 

4-5. Summary 

Let us summarize our results on the characteristic values of (-Bj) in the strip 
{z £ C : — 7r < < 7r} and introduce some notation for the following sections: 

• The set of characteristic values of (B%) in the strip {z £ C : — ir < SKz < it} 
is countable and symmetric with respect to the real axis. 

• An even number 2n of these characteristic values (counted with their to- 
tal null multiplicity) lies on the real axis. The corresponding Floquet 
modes are propagating, n linearly independent "outgoing" Floquet modes 
Vi,...v± are selected as described in section [2] or |Appcndix A| and the 
corresponding characteristic values (or quasi-momenta) are denoted by 

SI > • • • j S5 ' 

• We arrange the characteristic values with positive imaginary parts (counted 
with their total null multiplicities) in non-decreasing order of their imagi- 
nary parts and labeled them ^ +1 ,^ +2 , ■ ■ ■ ■ Corresponding exponentially 
decaying Floquet modes are chosen as in Proposition 13. 5p and denoted by 
id, v v£ +2 , .... For n > N the characteristic values are close to the 
reference values £^ and belong to the set S sketched in Fig. |2(b)[ 

• For Dirichlet and Neumann lateral boundary conditions the characteristic 
values for n > N are simple. For other boundary conditions the total 
null multiplicity of with n > JV is always < 2. and except for quasi- 
periodic boundary conditions with /3 £ {0, ir} at most a finite number of 
them is not simple. 
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5. Properties of the operator T 



5.1. Reference operator T T 

To study the operator T : (a n ) i— > Y]°°_ t a n v+ in (|13[) . we introduce a refer- 
ence operator T r using the Floquet modes 



\ -1/2 



v T n {x) := k, + — e-^^nC^), «£N if k„ > (23) 

corresponding to to — 0. The normalization constant has been chosen such that 
ll^nll h 1 (s+) = 1- in the special case k„ = 0, which can only occur for n = 1, we 
set Ui(x) := e _:ci -0i(a;2). In this case is not a solution to (fTaj) with w = 0, 
but we will not need this property. T r is defined by 

oo 

T r : Z 2 (N) -> i^< + (S+), (a n ) h> £ a„<, 

n=l 

and has the following properties: 

Lemma 5.1. 1. T r is weZi defined and isometric with T r (Z 2 (N)) C i?^(S l+ ) ; 
i.e. ||T r (a„)|| ff i = ||(o„)||p for all (a„) G l 2 (N). 
2. tT t : I (N) — > HZ is a Fredholm operator with index 0. 

Proof. Part 1: The assertion follows from the fact that {wj^ : n £ N} is an 
orthonormal system in H*' + (S + ). 

( i \ _1/2 

Fart I?: Define ip n := in n + ^— J V« if K n > and tp n := tp n else. Then 



(rT r )(a„) = X)^Li(^d — KnO^a-nipn- If #n = 0, the assertion follows from the 

r 

7 



~ 1 /2 

fact that {V'n : ?i 6 N} is a complete orthogonal system in HJ. = {(0,L)) 



and sup ngN H^nllfl-i/a/infneN HV'nllir 1 / 2 < 00 • To treat the case £?n note 

~ —1/2 

that {V'n : n G N} is also a complete orthogonal system in = H 1 ((0, L)) 

and that sup neN ((l + n 2 n ) 1/2 \\^ n \\ H -^^) I inf„ eN ((l + k^^HVVsIIh-vs) < oo. 
If none of the coefficients 6^ — K n 9^, n £ N vanishes, this immediately implies 
that t% : l 2 (N) -> Hy 1/2 ((0, L)) is bounded and boundedly invertible. If 
Ojy K n ^]sj = precisely for n = no, then the operator (a n ) i-> (TT r )(a„)+a rao '0 rl , o , 
which is a rank-1 perturbation of rT r , is boundedly invertible and hence tT v is 
Fredholm with index 0. The case that two of the coefficients 6*d — K n 9^ vanish 
can be treated analogously, and because of the form ([5]) of the n n no more than 
two coefficients can vanish. □ 

5.2. Estimates on perturbation of eigenvectors 

In the following estimates on the perturbation of eigenvectors we will make 
all constants explicit. To do so, we first have to introduce for each connected 
component S of S in (|2ip an L 2 -orthogonal projection operator 

P s u ■= J2 ^' ^o,n)Vo, n with I s := {n G N : £ G (24) 
nels 
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and define the quantity 

Sj := - inf{|«;„ — K n '\ ■ n,n' e N, re„ ^ K n >} ■ 

It is easy to see from the explicit values of the n n given in § l3.1l that <5 7 = ir/(2L) 
for 7 e {7d,7n,7dn}, S~, p = mhn eN \lir - /3|/i if /3 <£ {0,7r}, and 8 lft = n/L 
for /3 € {0, 7r}. Note that for a connected component 5 of the set 5 defined in 
Proposition 13.51 we have 

|9£ - 9£o,n| > ^7 for ah £ e 5, £ 5 (25) 

(see Fig. [3(b)]). 

Lemma 5.2 (estimates of eigenvectors). For a connected component S of S in 
(|2ip define Kg := inf^ g j 5£ and recall the notation (j24[) . TTien aZZ eigenvectors 
u* corresponding to characteristic values in S satisfy the estimates 

ifP'E 

- PsU*||z,2(n) < — t-tt r — |K||z,*(n), (26a) 

IU=(n) < lKllz=(n), (26b) 

2 - 

||(u' - P 5 u*)(0, OlliW)) < C^\\u*\\ L 2 (n) , (26c) 



<9u* 



9a; i 

withC:= £^-,1)1/3. 



9, ,2 = 

1, '' ) ,«,^ IWt » (2M) 



Proof. Let ( € 5 be a characteristic value of (Bj) with eigenvector u*. It follows 



from (HI that (£ - - f»)«*(0 = ^ 2 (e P w*)(0, or 



~ = ^ 2 (£ P u*)(/^ = , , , . 

This identity will be used extensively. Moreover, we need the lower bounds 

|e-6| 2 = (^ + 2-m) 2 + (^-.„) 2 >{[ 7rm JJ' A2 (28a) 

|£ - 6| 2 = («£ + 27T77i) 2 + (9£ + «n) 2 > (7T77i) 2 + (/SJ + K„) 2 (28b) 

which hold for all £ € 5 and (gZx N\ {0} x 7j and follow by separate estimation 
of real and imaginary parts using ([25]l . |9?£| < 7T, k„ > (see Fig. |3(b)j ). 
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Estimating — — > min(&y,7r)K£ forle^x N\ {0} x /§ we obtain 



lit* — 



5" Hi 2 



< 



min(<5 7 , n)tti 



E u*(i)(pi 

iGZxN\{0}x/5 
2 

E im*xoi 2 < 

zeZxN 



E 

;ezxN\{o}x/,5 



w 2 (e p M*)(0 



min(i5 7 , 7r)K^ 



Kill- 



To prove (|26b[) we use that |£ — £;| |£ — £;| > /cg^lm! and dxi^o,™) = to obtain 



|| L 2 



< 



E u*(l)d Xl <pi 

ZG(Z\{0})xN 
V «5 



= E 

/G(Z\{0})xN 



(27rm)w 2 (e p M*)(0 



I £3- 



Since 



we have 



s (o,o= E ^(o^(o,-) = E ( E u *( m > n )p»< 

ZGZxN n6N VmeZ / 



||(u*-PsU*)(0,-)||i a( ( Oli ))= E |E^( m ' n ) 

n£N\I s raeZ 
+ E J E M *("V 

ne/j mez\{o} 



(29) 



Using the Cauchy-Schwarz inequality, (f27|) . and the lower bounds |£ — £j| 2 > 
7r 2 m 2 + (5 2 and |£ - ^| 2 > k 2 ~ (see p5]l). the first term in ((29]) can be bounded 
by 



E 



E m *("v 



< 



( ^ 7r 2 m 2 + a ) ^ ^ 

VmeZ T/ mGZ,nGN\/j 



! m 2 +<5 2 )|M*(m,n)p 



<c 2 e 

meZ,neN\i,s 



w 4 i(e pM *)(™,n)i 2 <C 2 ^y n^ii 



4 



2 

L 2 • 



Using (|26b[) the second term in (f2"9")l can be estimated by 



E E u *( m > n ) 



< 



E — 

^ (27rm) 2 

vmeZ\{0} v ' 



E] (27r™) 2 |w*(m, n) 

v m<EZ\{0},rieJg 



C 2 



< T R 1 u*||i 2 <c- 



a (^ 2 „...„2 
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completing the proof of (|26c[) . 
To prove (|26d[) . first note that 



;ezxN neN \mez\{o} / 

It follows from (|27|l. the Cauchy-Schwarz inequality, the estimate |£— £/| |£— £j| > 
(■jrm) 2 (see ([2"5])). and the identity J2m=i m_2 = IT * na * 



\\d xl u*(0,-)\\h=J2 



iiei 



i 

2/ 



y-v L0 2 (2-Kim)(£ p u*){l) 



< £ (W f )2 ' ( /^ )|2 < 

\ ^-^ m z / ^— ' (7rm) 4 3 

\ m=l / (gZ\{0}xN v ' 

which proves (|26d|) . □ 

We can still fix a complex scaling constant in the Floquct modes v£ defined 
in §4.51 For our purposes it will be sufficient to fix this constant for all but a 
finite number of these w+. We will assume in the following that 

7^{7o,7j- (30) 

The special case 7 £ {70, 7?r} will be discussed in Appendix | Appendix C| If (|30|) 
holds true, it follows from the explicit form of n n given in § [3] that there exists 
N £ N such that for all n > N the connected component S n of S containing £Jj is 
a disk containing precisely one characteristic value with multiplicity. Therefore, 
the function 

Un(x) := e-^ xi v+(x), n>N (31) 

is periodic in x\, P§ u n is constant in the first variable, and (Pj u n )(0, ■) is a 
multiple of ip n . Therefore, we can choose the free complex scaling constants of 
v + and u n such that 

1 \ -V2 



(Ps n u„)(0, x 2 ) = U„ + — ?p n {x 2 ) . (32) 

If uj = this scaling yields = v T n . 

Corollary 5.3. There exist a constant C > depending only on tu 2 e and 7 
such that for all n > N 



*T K - P§ n un || £8(n) <C, (33a) 
"n-Ps n ^)(0,-)|| i2((0ii)) <C, (33b) 

K n 2 ll^x 1 WnUi=o|li>(( 0j i)) < (33c) 
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Proof. It follows from Lemma [5.21 that ||Mn||L 2 (Q) < 2||Pg Uti||l 2 (si) for suffi- 

ciently large n. Moreover, due to {52"} we have H-Pj^itrallz, 2 ^) = [ K n + ^T) 
Now the assertions follow from Lemma |572"1 □ 



5.3. Properties of the operator T 

Proposition 5.4 (properties of T). 1. The operator T is well-defined and 
bounded, and T — T r is compact. 
2. T(a n ) satisfies {la} and {lb} for all (a n ) G l 2 (N). 

Proof. Part 1: In the following C denotes a generic constant depending only on 
w 2 eand7. Note that v r n (x) = exp(— K ra xi)^> n (x2) with ipn(x2) '■= {Pg u n )(0,X2) 
for n > N. Inserting the term zte l ^ nXl ^j n (x2) in ||i>+ — i^lli 2 and applying the 
triangle inequality yields the estimate 

\K-<\\l'(S+) < \\^ X1 (u n -P Sn U n )\\ L2(s+) + \\e^-e-^\\ L2 ||VgU 2 



* + {J l^-^IJW- 

< -573, n > N. (34) 

In the second line we have used that < 3<^y < 3£„ for all n > N and 
J2T=o e" fca? " = (1- e"^")- 1 < (1 - e-'^f)- 1 , and in the third line {55a} was 
applied together with the identity 

I ri — k Ti 1^ 7 £n 

e u 1 - e ™ 1 dxi = 



-5£n(£n ^n) ^n(£n "i" i^n) 

and Proposition [473] 

To obtain an identity for the L 2 -distancc of the gradients, we apply Green's 
first theorem in (0, I) x (0, L), use the identity A (u+ — v r n ) — — uj 2 e p v+, and let 
I —} 00: 

||V(u+ - Vn)\\ 2 L 2 {s+) =w 2 J e P v+{v£ -v^jdx 

L 8(v+ -vl) X2) (-+_ ^r )(0j X2) dx2 

Since a ^"g~"" ) (0, ^2) = «6i u n(0, £2) + Kn^n + §|^(0, x 2 ), it follows after adding 
±i£,nipn and using {33b}, {33c}, and Proposition 14.31 that 
" ^LT^ MIU 2 - C/V^"- Usin 8' Cauchy's inequality, {34}, and {33b} yields 



0: 



V(tf+-i£)|| i2(S+) < — , n>iV. (35) 
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Define K s : V -> H}- + {S+) by := £ J n=1 a n {v+ - Combining (01 

and (|35|) and using Cauchy's inequality, we deduce that 

("12 , 7712 , \ 

E E ^ iimii 2 , 

which implies together with ([5]) that (A'j) is a Cauchy sequence with respect 
to the operator norm. Therefore, K = lim,,-^ Kj is well defined, and since 
the range of the operators Kj is finite dimensional, K is compact. Moreover, 
T = T r + K is well-defined and bounded. 

Part 2: Since the differential operator A + oj 2 e p is continuous from H*(S + ) 
to H~ 1 (S + ), we can interchange its application with summation to show that 
w := T(a n ) satisfies (|Ta|) . Analogously, it follows from the continuity of the 
trace operators that w satisfies (jlbl) . □ 



6. Proof of the Main Theorem 
Proposition 6.1. The operator T is infective. 

Proof. Let v := X^^Li a » w n f° r some sequence (a„) £ l 2 (N) and assume that 
v = 0. We have to show that a n = for all n <E N. Let {v\, v 2 , • ■ • } = : 
n e N} with < v\ < v% < ■ ■ ■ . We show by induction in I G N that a n = 
for all n € N satisfying 3£ n < v\: Formally adding v§ := —1, the induction 
base is trivial. For the induction step assume that the statement holds true for 
I — 1 and that there are precisely M distinct characteristic values £i, . . . £m with 
S£ m = v i an d £ [— 7r, tt)- Pccularities for the case I = 1 and vi = will 
be discussed at the end of the proof, so assume for the moment that v\ > 0. 
Suppose that n((i?{), £ m ) = (r^ , . . . , ra^ ). Relabelling the Floquet modes i>+ 
and the coefficients a n corresponding to £i, . . . , by and we obtain 
from the induction hypothesis that for any e > and all x 2 £ [0, L] 

exp^r)^^ a $ v $( x ) = °( e ~ i ' n+1 ~ vi ~ e)xi )> xi ^oo (36) 

m=l j'=l fc=0 

and have to show that all coefficients vanish. For simplicity, let us first as- 
sume that all partial null multiplicities r^ are equal to f . With the notation of 
Proposition ^. 51 fand an additional index m) the function exp(— i£ m %i) 5Z"=i a< fo ^ v< jo' \ x ) — 

a jfl^o}n( x ) * s 1-periodic in x\ for each m. Pick some x = {x\,X2) € H 
and set g m := X)j=i "p^ait 3 ")' Then it follows from (f36|) that 

M M a m 

^ cxp(^| m ) p 5m = exp(^p) ^ exp(«| ro p) ^ affiu^ m (x) (37) 

m— 1 m— 1 J — 1 

1/ Q m 

= exp(i^p) E ?'™^!^! +P' :E 2) -> 0, as p -> oo ,p e N. 

ro=l j=l 
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If the right hand side would vanish exactly for M consecutive values of p, say 
p £ {<?+l, ■ • • , q+M}, we could conclude immediately that g := (gi, . . . , <?a/) T G 
C M is zero since the matrix e C MxA/ defined by := exp(i^ m ) / + < ? ! 
I.m = 1, ...,M is regular as shown below. However, (|37[) only implies that 
lim g _ i . 00 ^4( 9 )g = 0. Therefore, we have to control || L4( ? )] -1 1| uniformly in q. For 
this end note that A^ has a factorization 

A (i) = A (o) diag ^(i^!) 9 , . . . , exp^A/) 9 ) 

and ^4 ( - ^ is a Vandermonde matrix. It follows that A^ is invertible and 
|| ^("J)] -1 ] is independent of q. Therefore, lim g _ i . 00 A^g = implies g = 0, 
i.e. a j o m = for each m = 1, . . . , M . Since we know from Proposition 

13.51 that UQ m , . . . jWom^ arG linearly independent, it follows that a\ "q = ■ ■ ■ = 

Now assume that f := max{rj m ' ) : m = 1, . . . , M,j = 1, . . . , a^} is greater 
than 1. Multiplying (J36J) by xT/ 4-1 and using ([II]) it follows that ([37]) holds 
with <7 TO := X)j. r < m )-r a j^-i"o'm( a ') ( w hh the convention that an empty sum is 
0), and it follows as above that g\ = ■ ■ ■ = g M = for all x £ fl. Using again 
the linear indepence of the functions v^ m , we find that a) > ™_ 1 = for all (j, m) 
such that = r. In a second step we show analogously that o^_ 2 = ^ or 

all (j,m) such that > r — 1, and so on. Finally, all coefficients ay^ have 
to vanish. 

It remains to discuss the case of propagating modes, i.e. the induction step 
for z^i = since not all elements of the eigenspaces are considered here. How- 
ever, we can use the same technique as above to prove the stronger result that 
En=i a n v n +Y,n=i a t v i = implies = a{ = ■ ■ ■ = a- = a+ = a+ = ■ ■ • . □ 

The well posedness result in the following proposition follows from general 
results in [19j . but we obtain an independent proof as a side product of our 
analysis: 

Proposition 6.2 (well-posedncss of problem ([I])). Assume that the only solu- 
tion v £ H} f ,+ (S + ) to |I]) with f = is V = 0. Then F := tT is bounded 
and boundedly invertible, and problem |T]) is well posed in the sense that for 
all f £ there exists a unique solution v £ H} / ,+ (S + ) to ([T]), and v depends 
continuously on f. 

Proof. F is a Fredholm operator with index since F = tT v + t(T — T r ), and rT r 
is Fredholm with index by Lemma l5.ll and t(T— T r ) is compact by Proposition 
15.41 Assume that F(a n ) = and set v := T(a n ). Then v is a solution to ([TJ 
with / = 0, and hence v = by our assumption. Using Proposition 16.11 we 
conclude that (a n ) = 0, i.e. F is injective. This implies that F has a bounded 
inverse. Using Proposition 15 Al it follows that v := TF~ 1 f is a solution of ([T]), 
which depends continuously on /. □ 
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The proof of our main theorem is now simple: 

Proof of Theorem ] 2. SI We start with part 2 of the theorem: Since F := tT is 
bounded and boundcdly invertiblc by Proposition l6.2[ the set {ru+ : n £ N} is a 
Riesz basis of H'L. It follows from Propositions l6 . 21 and IB .21 that ([T]) is well posed 
for tv = v(0, ■) + z^j-(0, •). For this choice of r, the operator T has the bounded 
left inverse F _1 t, and hence {v+ : n £ N} is a Riesz basis of ran(T) = V. 

Since the functions v„ are chosen as in Proposition 13.51 the matrix repre- 
senting T consists of Jordan blocks. Because of Proposition 14.31 at most a finite 
number of these Jordan blocks have size > 1. It is straightforward to see that 
TZ is represented by the same matrix. □ 



Appendix A. On the radiation condition 

The study of solutions to the Hclmholtz equation Av + u> 2 e p v = in S with 
boundary conditions jv(xi,-) = (0, 0) T for x\ £ M amounts to the study of 
spectral properties of the operator := ——A : H 2 (S) — s- L 2 (S). Due to 
the isometry of the Floquet transform, the spectrum of A 1 is the union of the 
spectra of the operators defined by 

:= - — A a : ff^(O) -> L 2 (n), for «el. 

Since these operators are positive and self-adjoint in the weighted Hilbert space 
L 2 (f2, e p ) with a compact resolvent, their spectra consist of a countable number 
of positive eigenvalues with finite multiplicities accumulating only at oo: 

a(A^) = {\ m (a) : m £ N], a£ [-tt, tt). 

We assume that the A m arc arranged in increasing order. The functions A m 
are smooth, except at points where two or more of them cross. Alternatively, 
the eigenvalues of can be arranged such that they are analytic functions 
A rn of a, which have holomorphic extensions to a complex neighborhood U of 
[— 7r, tt) ((3). Furthermore, there exists a holomorphic family of eigenfunctions 
U L 2 (il), £ i-> w n ^: 

A^w m ^ + A m (0e P Wm,€ =0, (6«,meN. (A.l) 
Note that if £* £ [—tt, tt) is a characteristic value of (Bg), then 

kcrB c » = ker(A^ } - w 2 I) = span{w m ^» : 3m £ N X m (C) = uj2 }- ( A - 2 ) 

Proposition A.l. £* £ [— tt,tt) is a characteristic value of (B^) with a partial 
null multiplicity greater than 1 if and only if there exists m £ N such that 
\ m (C)=0J 2 and\' m (e)=Q- 
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Proof. Suppose that uj 2 = A m (£*) and X' m {C) = 0. Taking the derivative of 
(|A.1[) with respect to £, which will be indicated by a prime in the rest of this 
proof, yields 



(A e + A m (£)e P ) w'e + 2i(d xi + iQw m ,t + A^(£)e P w m , e = 0. 



(A.3) 



Since \' m {£,*) = 0, we have 



o 



and 



BfWm^* = 



Therefore, is a root function of (_Bt) corresponding to £* with a partial 

null multiplicity greater than 1. 

Conversely, assume that £* is a characteristic value of B% with a partial null 
multiplicity greater than 1. Then, there exists a root function (itj) such that 
B{*U£* = and 

B^u't. +B' i ,u^=0. (A.4) 

Let S(cj 2 ,£*) := {m e N : A m (£*) = w 2 }. Due to (|A.2[) there exist coefficients 
v rn e C such that = ug» with u>£ := X)mes(w 2 5*) v mW m ,i- Taking a linear 
combination of the equations (|A.3[) and subtracting (|A.4[) . we get 



E A ™(^) 

m£EH(a; 2 ,£;*) 



The right hand side of (|A.5[) belongs to ran(i?£«). Since £* e 
adjoint in L 2 (f2), and hence ker(B^*) = ran(B^»)- L . As 

E Kn(C)^mW m ,^ € ker(.B f .)> 



(A.5) 
i?^* is self- 



we have 



E A »(n 



= 0. 



As the functions w m t* are linearly independent, it follows that A^ (£*);/,„ = 
for all m £ S(w 2 ,£*). Since not all v m vanish, we obtain that A^(£*) = for 
some m. □ 



Recall that the group velocity of a Floquet mode w m ^(x)e l ^ Xl is given by 

(A.6) 



di 



Proposition A. 2. For a Floquet mode of the form v(x) = W m ,£*{x) 
e 1 ^ X1 , A TO (£*) ~ with non vanishing group velocity, i.e. X' m (£*) 7^ the 
following statements are equivalent: 
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1. v has positive energy flux, i.e. ^sq(v,v) > 0. 

2. v has positive group velocity, i.e. \' m {£,*) > 0. 

Moreover, if v(x) = w n ^*(x)e 1 ^ Xl is another Floquet mode with A„(£*) = uj 2 
and m ^ n, then q(v, v) = 0. 

Proof. Since q Xl (v,v) = q(v,v) is independent of x\, it follows that 
ni r dv 



q{v,v) 



dv 

dx ^x)v(x)-v(x) — (x) 



dxo 



(x)v(x) - v(x)- — (x) 



o Jo 



dx\ 



dx\ 



dx2 dx\ 



dx. 



Taking the L 2 inner product with w n .£* in (|A.3[) and using (B^w' m £,,w n ,$*) = 
(w' m £», B^*w n ^) = and analogously with the roles of w m ^* and u>„,£* inter- 
changed we obtain that 



' dx. 



Since w m ^ and w n ^ arc orthogonal with respect to the inner product in 
L 2 (f2,£ p ) for all £ G K as cigenfunctions of the self-adjoint operators ^ in 
this space, we obtain the last statement. Choosing v = v shows the equivalence 
result. □ 

Remark A. 3. Two further important items could be added to the list of equiv- 
alent statements in Proposition IA.2 I which we do not want to define in detail 
here: It has been shown by Fliss [ly, Theorem 3.2.57] that a Floquet mode 
has positive group velocity if and only if it satisfies the limit absorption prin- 
ciple. Moreover, the equivalence of the principles of limit absorption and limit 
amplitude has been shown in the thesis of Radosz [23} . These are very strong 
indications that solutions with positive group velocity are "physical solutions" , 
and Proposition IA. 21 shows that for frequencies uj for which all real characteris- 
tic values have total multiplicity 1 (cf. Proposition lA.il for a characterization), 
the conditions ([3]) are satisfied precisely for the physical Floquet modes (up to 
normalization) . 

However, if two (or more) bands X m and A„ cross at oj 2 , i.e. A m (£*) = 
A n (£*) = uj 2 for some £* e [— 7r,7r) and if AJ 71 (£*)A^ 1 (£*) < 0, then a system of 
Floquet modes satisfying @ docs not necessarily satisfy the limiting absorption 
principle. E.g., the Floquet modes v a :— V2v + v, Vb := y/2v + v satisfy (j3)), but 
v a does not satisfy the limiting absorption principle. This shows in particular 
that Conjecture 4.2 in [24[ is false in general, but true if all real characteristic 
values have total multiplicity 1. 

For Floquet modes with group velocity we do not know which are the 
physical solutions. A system of Floquet modes satisfying (|3|) is constructed 
in (19( 1 , but we are not aware of indications that they correspond to physical 
solutions. 
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Appendix B. Uniqueness results 

In this appendix we discuss conditions under which solutions to the boundary 
value problem ([TJ are unique. In this case problem ([1} is well posed (Proposition 
16.21) . and the second part of our main theorem 12.21 holds true. Proposition IB. 21 
is used in the proof of Theorem 12.21 

Lemma B.l. Assume that either 8^ = (Dirichlet condition) or (9n = 1 
and > 0. Moreover, assume that the only solution v G -ff*(S ,+ ) to the 
boundary value problem (jla|) - (jlcp with f = is v = 0. Then the only solution 
v G H^ + {S + ) to (Ta}-([Tc]) with f = 0isv = 0. 

Proof. Let v = J2n=i CnV n + w £ -^7' + ( , S' + ) De a solution to (jTa])- (JTcJ) with 
/ = 0. Then it follows from (lil . Theorem 5.1.4] that there exists 6 > such 
that exp(fei)ui G H} f (S + ). Therefore, letting x\ tend to oo in the definition of 
q(w,v^), it follows that q(w,v^) = for all j, and analogously q(w,w) = 0. 
For the case #n = 1 and > Owe obtain from ^-(0,^2) = -<?d"(0,^2) and 
the orthogonality and normalization conditions for the Floquet modes v£ that 

O>-2(30 D )/ K0,a; 2 )| 2 da; 2 = 23 / — (0, :r 2 )y(0, x 2 ) cfo 2 = 
Jo Jo "^l 

= EM 2 - 

3=1 

This implies c\ = ■ ■ ■ = c— = and hence v = w. For the case #n = it can 
be shown analogously that v = w. Now the assumption of the lemma implies 
v = 0. □ 

Proposition B.2. Let tv = tJ^(0,-) + inv(0, ■) with n > 0. TTien i/ie only 
solution v G to (P wit/i / = is v = 0. 

Proof. According to Lemma IB . 1 1 it suffices to consider solutions v to ([I} with 
/ = in H} ( (S + ). Then there exists a sequence tending to 00 such that 

lirn^oo I J L -§^{Rk, x 2 )v(R k ,x 2 ) dx 2 \ = 0, and hence q(v, v) = 0. Therefore, 

f L dv f L 
= $q(v, v) = 29 / —(0, xsMO, x 2 ) dx 2 = 2k / \v(0, x 2 )\ 2 dx 2 . 
Jo 0x1 J 

Hence, i>(0, •) = 0, and ^-(0, •) = zkv(0, •) = 0. Now a unique continuation 
principle in two dimensions (see [2^, Corollary 7.4.2]) implies that v = 0. □ 

Proposition B.3. Assume that t is either the Dirichlet or the Neumann trace 
and that £ p satisfies the symmetry condition e p (l — xi,x 2 ) = e p (xi,x 2 ). Then 
the only solution v G H} f ' + (S + ) to (fTJ) with f = is v = 0. 
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Proof. In both cases Lemma TB.ll applies, so it suffices to show that a solution 
v G H} ( (S + ) to flTJ with / = vanishes. For the Dirichlet condition u(0, •) = 
(i.e. 9b = 1 and (9n = 0) we consider the odd extension v(— x\, X2) '■= —v(xx, X2) 
for x\ > 0, X2 G (0, i), and for the Neumann condition J^-(0, ■) = (i.e. # N = 1 
and On — 0) the even extension v(—xi,X2) ■= v(xi,X2). In both cases the 
extended function v satisfies the differential equation (fTaj) in the full strip S, 
and hence B a (Fv)(-, a) = for all a. Since (B^) has at most a finite number of 
characteristic values on the real axis, it follows that (Fv)(-,a) = for almost 
all a. As v G L 2 (S), an application of the inverse Floquet transform yields 
v = 0. □ 



Appendix C. Quasi-periodic boundary conditions with (3 E {0, tv} 

If 7 G {70, 77r}, there are infinitely many connected components of S contain- 
ing two characteristic values (with multiplicities). Here we set N := N. Then 
each characteristic value with n > N shares a connected component S n of S 
with precisely one other characteristic value, w.l.o.g. For 7 G {70,7^} wc 

cannot exclude the possibility of infinitely many characteristic values with par- 
tial null multiplicity 2 in general. Treating this case would require an extension 
of Lemma 15.21 and considerable additional work. To get along with Lemma 15.21 
here, we have imposed the additional assumption e v (xi,X2) = e p (xi,L — X2) 
in the case 7 G {70,771-} above. Then by a symmetry argument we can split 
the problem into two problems for wave guides of width L/2: For /3 = wc 
either impose Dirichlet conditions at {x : X2 G {0, L/2, L}} or Neumann condi- 
tions at {x : X2 G {0, L/2,L}}, and for /? = tt we impose Dirichlet conditions 
at X2 = 0, X2 = L and Neumann conditions at X2 = L/2, or Neumann con- 
ditions at X2 = 0, X2 = L and Dirichlet conditions at X2 = L/2. For each of 
the subproblems for waveguides of width L/2 wc can exclude that possibility 
of infinitely many characteristic values with partial null multiplicities > 2 as 
above. Therefore, u n defined by pip is again periodic in x\. However, P§ u n 
is not necessarily a multiple of tj) n , but we only have P§ u n G span{-0„, ^> n ±i}. 
Therefore, we replace (|32[) by the normalization condition 

/ l \-Va 

\\Ps n u nh*{n) = ( K « + 7^r J . n>N, (C.l) 

which leaves a free phase factor. Moreover, we replace (f2Uf by 

v c n (x):=e- K ^{P Sn u n )(x), n>N. 

It is easy to see that Lemma I5TT1 and Corollary 15.31 remain valid. The proofs of 
other results are not affected. 
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